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Abstract 

In most models of traffic flow, the car density p is the only free parameter in 
determining the average car velocity (v). The critical car density p c , which is defined 
to be the car density separating the jamming phase (with (v) = 0) and the moving 
phase (with (v) > 0), is an important physical quantity to investigate. By means of 
simple statistical argument, we show that p c < 1 for the Biham-Middleton-Levine 
model of traffic flow in two or higher spatial dimensions. In particular, we show 
that p c < 11/12 in 2 dimension and p c <l - i^rr) in -D (D > 2) dimensions. 
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Recently, there has been much interest in the study of traffic flow problems within the context 
of cellular automaton (CA) models. This approach has the advantages of being simple and compu- 
tationally less intensive. Nagel and Schreckenberg |]J introduced a stochastic discrete automaton 
model in one dimension (ID) to study the transition from laminar traffic flow to start-to-stop-waves 
as the car density increases. Variations of the ID model have been extensively studied Biham, 
Middleton, and Levine (henceforth referred to as BML) introduced a simple two dimensional 
(2D) CA model with traffic lights and studied the average velocity of CclXS clS db function of their 
concentration. In the BML model, each site on a 2D square lattice may take on an east-bound 
car (i.e., a car traveling in positive x-direction) and a north-bound car (i.e., a car traveling in pos- 
itive y-direction) both with probability p/2. Otherwise, the site is empty with probability 1 — p. 
The probability p thus corresponds to the car density. Each east-bound car moves at even time 
steps by one site provided that the neighboring site to the east is empty. Similarly, each north- 
bound car moves at odd time steps by one site provided that the neighboring site to the north is 
empty. BML found that, using numerical simulations, the average car velocity in the long time 
limit vanishes when the car density is higher than a critical value p c . The BML model has been 
generalized [0-0, for example, to include effects of two-level crossing (overpasses), faulty traffic 
lights, accidents, and anisotropy in car densities in east and north bound directions. Generalization 
of the model to other spatial dimensions is straight forward. "Traffic" models in higher dimensions, 
though may not be relevant to cars on the road, may be related to problems in information network 
and telecommunication. 

Nagatani f?|,|8|] has proposed a mean field theory for the average car velocity. The theory has 
recently been criticized by Ishibashi and Fukui [11]]. They suggested that no complete traffic jam 
occurs unless p — 1, in contrast to simulation results and other mean field approximations. They 
also speculated that the values of p c obtained in numerical experiments may be results of finite size 



effects and the use of periodic boundary conditions [0. It is therefore desirable to derive bounds 



for p c in various spatial dimensions. In particular, if the upper bound of p c is less than unity, then 
there exists a range of car densities corresponding to a complete traffic jam. In what follows, we 
argue that there exists an upper bound for p c less than unity whenever the spatial dimension of the 
lattice is greater than one. In 2D, we have p c < 11/12 and we also give upper bounds for p c for 
spatial dimensions D > 2. 

Consider the BML traffic flow problem on a 2D square lattice. Instead of looking at how the 
traffic flows, we concentrate on how the traffic jams. We note that the 2D square lattice can be 
partitioned, for example, by a collection of inverted-L-shaped elements as shown in Fig. 1. Each 
element consists of 3 sites. For example, the element at the site contains a collection of three 
sites Lij = {i — 1, j), — 1)}. Thus, we can partition the 2D lattice by {L 3i+ jj : i,j 6 Z} 

(see Fig.l). Obviously, this is not the only way to partition the lattice. As we shall see, different 
partitions may lead to different upper bounds. It is, however, suffice for our purpose to consider 
this simple partition. The important point is that the elements can fill the underlying lattice with 
a density of non-zero measure. 

Consider the case of a full lattice of cars, i.e., p — 1, which trivially corresponds to a jamming 
phase. Consider, say, the element at (i, j). The car at site can be removed without changing the 
system to a moving phase if the cars at sites — and — 1) are north-bound and east-bound 
respectively. Assume the cars to have equal probability of being east or north bound, the probability 
of having the car at the middle of the inverted-L-shaped element removed without changing the 
system to a moving phase is (1/2) 2 . Since each of the inverted-L-shaped element contains three 
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sites, so in the limit of infinite system size, we can almost surely remove (1/3) x (1/2) 2 = 1/12 of 
the cars in the system for any initial configuration of a full lattice of cars without changing the 
system from a jamming phase to a moving phase. This suggests an upper bound of p c < 11/12 in 
2D. Conversely, let us consider the case for an arbitrary initial configuration with car concentration 
p — 11/12. A car moves into an empty site (i, j) in the next time step or two whenever there is a 
north-bound car located at — or an east-bound car located at — at the present moment. 
The probability for this to occur is given by g = y| | + if \ ~ (if f) 2 = I§ x H" Therefore, in 
the mean field approximation, the expected number of sites at a given empty site can move is 
given by 

(n)«(l-g)xO + g(l-g)xl + g 2 (l - q) x 2 H h q k (l - q) x k H 

So on the average, an empty site will be "blocked" after moving over (n) ~ 407/169 ~ 2.4 sites. This 
is less than the typical separation between two empty sites (which equals \f\2 « 3.4). Therefore, it 
is almost certain that a system with car concentration p = 11/12 will evolve eventually to a jamming 



state. Our result should be contrasted with that of ||11|| , in which it was found that complete jam 
occurs only when p — 1. We believe that the discrepancy is due to the fact that Ishibashi and Fukui 
have incorrectly assumed that all car configurations are independent of each other. The upper 
bound on p c that we obtained here is consistent with all the numerical experiments (Biham et al. 
found that p c ~ 0.32 for 512 x 512 lattice 0, while a finite size scaling analysis by Nagatani showed 
that p c = 0.42 J7|). We also note that our argument above is true for any boundary conditions. 



And a more careful choice of partitions and elimination procedure may lead to a tighter bound [|T2 

Our argument can be readily generalized to anisotropic distribution of north and east bound 
cars as long as the traffic problem is not decoupled into independent ID lanes. In ID, it is known 
that complete traffic jam only occurs when p = 1 Q. Let q e (q n ) be the probability of finding an 
east-bound (north-bound) car at a site given that there is a car at the site. We have q e + q n = 1, 
and the ratio q e /q n reflects the anisotropy. In this case, our argument leads to a bound in 2D 

qeqn fn \ 

Pc < 1 - -jp (2) 

Even in the ID limit, either q e or q n vanishes and p c < 1, which is consistent with the known results 
Pc = 1. 

Our argument can be carried over to any spatial dimension D > 2. A naive extension of the 
above partition to 3D may be that of introducing elements with four sites. For example, the element 
Lij k at site (i, j, k) contains the collection of sites L^. = k), — k), — k), k — 1)}. 
However, one cannot partition a 3D cubic lattice with elements Since the choice of partition 
is quite arbitrary, we choose, for simplicity, to fill the underlying lattice with hypercubic elements 
in D dimension, each of which contains 2 D sites. Consider a full lattice of cars. The probability of 
removing a car from a hypercube without changing the jamming character is (1 — j=j) D '. Thus for 
D > 2, we have an upper bound for p c given by 

p ^ 1 "^( 1 -i)) D = 1 -(^ir) D - (3 > 

In 3D, this gives p c < 26/27. For anisotropic distributions in different directions, we have 
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Pc <l-^f[(l-q e ), (4) 

t=\ 

where qe is the probability of finding 8b CclX ell db Sit e moving in the £-th direction provided that there 
is a car at the site, and thus qe — 1. As a result, in the infinite dimensional limit, we find 
that p c < 1 which is consistent with the mean field theory result that (v) — 1 — p where p is the 
probability of finding a car in a site. 

In summary, we argued that the critical car density has an upper bound less then unity in traffic 
flow problems in D > 2. 
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FIGURES 

FIG. 1. Partition of a two dimensional square lattice by elements containing three sites. The dashed 
line outline one of the elements. 
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Fig. 1 



